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Dedicated with great pleasure to Ludwig Streit on the occasion of his 60th birthday 

Abstract. An elementary yet remarkable similarity between the Cole-Hopf 
transformation relating the Burgers and heat equation and Miura's transfor- 
mation connecting the KdV and mKdV equations is studied in detail. 



1. Introduction 

Our aim in this note is to display the close similarity between the well-known 
Cole-Hopf transformation relating the Burgers and the heat equation, and the cele- 
brated Miura transform connecting the Korteweg-de Vries (KdV) and the modified 
KdV (mKdV) equation. In doi ng so we wi ll introduce an additional twist in the 
Cole-Hopf transformation (cf. ( 1.28 ), ( 1.2S| )), which to the best of our knowledge, 



appears to be new. Moreover, we will reveal the history of this transformation and 
uncover several instances of its rediscovery (including those by Cole and Hopf). 

We start with a brief introductory account on the KdV and mKdV equations. 
The KdV equation p3] was derived as an equation modeling the behavior of shal- 
low water waves moving in one direction by Korteweg and his student de Vries 
in 1895Q. The landmark discovery of the inverse scattering method by Gardner, 
Green, Kruskal, and Miura in 1967 |ll| (cf. also ^0|) brought the KdV equation 
to the forefront of mathematical physics, and started the phenomenal development 
involving multiple disciplines of science as well as several branches of mathematics. 

The KdV equation (in a setting convenient for our purpose) reads 

KdV(^) =V t - 6VV X + V xxx = 0, (1.1) 
while its modified counterpart, the mKdV equation, equals 

mKdV(0) = cf> t - &4> 2 4>x + 4>xxx = 0. (1.2) 
Miura's fundamental discovery Eq] was the realization that if 4> satisfies the mKdV 



equation ( |l.2| ), then 

V±(x,t) = <f>{x,t) 2 ±<j> x (x,t), (ii)el 2 (1.3) 



both satisfy the KdV equation. The transformation (1.3) has since been called 
the Miura transformation. Furthermore, explicit calculations by Miura showed the 
validity of the identity 

KdV(V±) = (20 ± d x ) mKdV(^). (1.4) 



Supported in part by the Research Council of Norway under grant 107510/410 and the Uni- 
versity of Missouri Research Board grant RB-97-086. 

1 But the equation had been derived earlier by Boussinesq [Q in 1871, see Heyerhoff and 
Pego [g. 
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The Miura transformation ( |l.3| ) was quite prominently used in the construction of 
an infinite series of conservation laws for the KdV equation, see [^J, Sect. 5.1. 
Miura's identity ( fOl ) then demonstrates how to transfer solutions of the mKdV 
equation to solutions of the KdV equation, but due to the nontrivial kernel of 
(2(j) ± <9 X ), it is not immediately clear how to reverse the procedure and to transfer 
solutions of the KdV equation to solutions of the mKdV equation. It was shown in 
(24) (sec also Q, [Q, |(|, |l|) how to revert the process. 

Following a similar treatment of the (modified) Kadomtsev-Petviashvili equation 
in 27|], we use here a method that considerably simplies the proofs in [^l| , , |^6| , 



|31 . Introduce the first-order differential expression 

P(V) = 2V8 X -V X . (1.5) 

Then one derives 

mKdV(0) = d x f ityt - PO4)V0j , (1-6) 

where 

^ = ^M ^>o, v± = 4> 2 ±4> x . (1.7) 

Next, let V = V(x, t) be a solution of the KdV equation, KdV(V) = 0, and ijj > 
be a function satisfying 

i> t =P(V)iP, -i> xx + ViP = 0. (1.8) 

Then one immediately deduces that c/> solves the mKdV equation, mKdV(</>) = 0, 
and hence the Miura transformation has been "inverted" . 



The KdV equation (LI) and the mKdV equation (L2) are just the first (non- 



linear) evolution equations in a countably infinite hierarchy of such equations (the 



(m)KdV hierarchy). The considerations (1.3)-( |1.8| ) extend to the entire hierarchy 



of these equations, replacing the first-order differential expression P(V) — Pi(V) 
by an appropriate first-order differential expression P n (V) for n G N (cf., e.g., 



21 1, p§, M, pi). More precisely, denoting the nth KdV equation in the KdV 



hierarchy by 

KdV„(V) = 0, (1.9) 

Lax (i^] constructed differential expressions P2 n +i(V) of order 2n + 1 with coeffi- 
cients differential polynomials of V such that 

j t L-[P 2n+1 {V),L] =KdV„00, neN. (1.10) 
Here L denotes the Schrodinger differential expression 

L = -d 2 x + V. (1.11) 



The KdV functional in (1.1) then corresponds to n = 1 and one obtains 

P 3 (V) = -4dl + 6Vd x + 3V x (1.12) 

in this case. Restriction of P2n+i{V) to the (algebraic) nullspace of L then yields 
the first-order differential expression 

P n (V) = P 2 n + i(V)\ kci{Ly n G N. (1.13) 
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Next we turn to the the Cole-Hopf transformation and its history. The classical 
Cole-Hopf transformation |l3| , |4lJ , covered in most textbooks on partial differen- 
tial equations, states that 

V(x,t) = -2^^-, (M)eRx(0,oo), (1.14) 

where tp > is a solution of the heat equation 

ipt=4>xx, (1-15) 

satisfies the (viscous) Burgers equation 

V t + VV x = V xx . (1.16) 

However, already in 1906, Forsyth, in his multi- volume treatise on differential equa- 
tions ( , p. 100), discussed the equation (in his notation) 

d fl / da , \ ] n da „ ,„ „. 

d- x Uv^- a )\-% = ^ (L17) 

where a = a(x,y). Hence there exists a function 9 such that 

90 da 2 00 ,„ oN 

a= 7T' 7-7r-a 2 = 2/3— . 1.18 
ox ox oy 

Assuming the function z satisfies 

z xx + 2az x + 2f3z y + ~f = Q, (1.19) 
an easy calculation shows that 

^ {ze o ) + 2 ^ (ze e ) = Q _ (1 . 20) 

Introducing new variables t = — y and u(x, t) = — 2a(x, y) as well as fixing f3 = 1/2, 
7 = 0, and z = 1, one concludes that (1.17) indeed reduces to the viscous Burgers 
equation 

u t + uu x = u XXl (1-21) 

while ( 1.2C| ) equals 

(e e )t = (e e ) xx , (1.22) 

with solutions related by 

u=-20 x . (1.23) 

However, Forsyth did not study the ramifications of this transformation, and no 
applications are discussed. 

Shortly thereafter, in 1915, Bateman Q introduced the model equation 

u t + uu x = vu xx . (1-24) 

He was interested in the vanishing viscosity limit, that is, when v — > 0. By studying 
solutions of the form u = F(x + Ut), he concluded that "the question of the limiting 
form of the motion of a viscous fluid when the viscosity tends to zero requires very 
careful investigation" . 
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Only in 1940 did Burgers (Q, p. 8) introduce^] what has later been called the 
(viscous) Burgers equation, as a simple model of turbulence, and did some prelim- 
inary investigation on properties of the solution. 

Taking advantage of the later rediscovered Forsyth transformation by Cole and 
Hopf, Burgers continued the investigations of what he called the nonlinear diffusion 
equation, focusing mainly on statistical aspects of the equation. The results of 
these investigations were collected in his book |ll| . 

In 1948, Florin |l7j], in the context of applications to watersaturated flow, redis- 
covered Forsyth's transformation, which would become well-known under the name 
Cole-Hopf transformation only some 44 years later. 

Although the Cole-Hopf transformation had already been published in 1906, it 
was only with the seminal papers by Hopf pl|P| in 1950Q and by Cole (ll| in 1951 
that the full impact of the simple transformation was seen. In particular the careful 
study by Hopf concerning the vanishing viscosity limit represented a landmark in 
the emerging theory of conservation laws. Although the Cole-Hopf transformation 
is restricted to the Burgers equation, the insight and the motivation from this 
analysis has been of fundamental importance in the theory of conservation laws. 
Furthermore, Cole states the generalization of the Cole-Hopf transformation to 
a particular multi-dimensional system. More precisely, if tj) = ?p(x,t), (x,i) £ 
K™ x (0, oo), satisfies the n-dimensional heat equation 

ipt = vAip, v > 0, (1.25) 

and one defines 

V = -22/Vln(V0, (1-26) 

then V satisfies 

V t + (V ■ V)V = vAV, (1.27) 

and the vector-valued function V = V(x,t) £ 1" has as many components (i.e., 
n) as the dimension of the underlying space. Observe, in particular, that V is 
irrotational (i.e., V = VW for some W, or equivalently, curll^ = 0,). The multi- 
dimensional extension was rediscovered by Kuznetsov and Rozhdestvenskii p3j in 
1961. 

In this note we show the following relations, 

-d + — 

ip x tp 2 



V t + VV X - vV xx = 2v ( --d x + Z|) - y^ xx ) (1.28) 



= -2vd x (±-(ip t - vip xx )j , (1.29) 
whenever V = —2vip x /il) for a positive function ip. This clearly displays the nature 



of the Cole-Hopf transformation and closely resembles Miura's identity (1.4) and 
the relation (1.6). Even though identities (1.28) and ( 1.29] ) are elementary obser- 



vations, much to our surprise, they appear to have escaped notice in the extensive 



2 Frequently Burgers equations is quoted from his 1948 paper po[, but he had already intro- 
duced it in 1940. 

3 With a misprint in the title, writing Ut + uu x = fi xx rather than ut + uu x = \iu xx . 

4 Hopf |4l[ states in a footnote (p. 202) that he had the "Cole-Hopf transformation" already in 
1946, but "it was not until 1949 that I became sufficiently acq uaint ed with the recent development 
of fluid dynamics to be convinced that a theory based on ( |l.2^ ) could serve as an instructive 
introduction into some of the mathematical problems involved" . 
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literature on the Cole-Hopf transformation thus far. While both the KdV and 
mKdV equations are nonlinear partial differential equations, the case of the Burg- 
ers and heat equations just considered is a bit different since it relates a nonlinear 
and a linear partial differential equation (see also ||, Sect. 6.4). 

One can also extend the Cole-Hopf transformation to the case of a potential 
term F in the heat equation, see, for instance, |57|. Here the relation ( 1.29| ) reads 
as follows, 

V t + VV X - vV xx + 2vF x = -2vd x - vijj xx - Fi/j) \ , (1.30) 

whenever V = —2vd x In^ip) for a positive function tp. The case of Burgers' equation 
externally driven by a random potential term recently generated particular interest, 
see, for instance, [gj, ||], |55), p7| , |39| and the references therein. We also 
mention a very interesting application of the Cole-Hopf transformation to the pair 
of the telegraph and a nonlinear Boltzmann equation in |fi"o| , generalizing the pair 
of the heat and Burgers equation considered in this note. 



Equation ( 1.30 ) extends to the multi-dimensional case corresponding to ( 1.27 ) 



and one obtains 

V t +a(V- V)V- vAV + —VF = -— V ( -U t - vAi\) - Fip)) , (1.31) 

a a \yj v '/ 

whenever a G R\{0} and V = —(2i//a)\7\n(ip) for a positive function ip. 

Obviously there is a close similarity between the heat and the Burgers equation 



expressed by (1.21s), and Miura's identity (1.4) relating the mKdV and the KdV 
equation. The principal idea underlying these considerations being that one (hierar- 
chy of) evolutions equation(s) can be represented as a linear differential expression 
acting on another (hierarchy of) evolution equation(s). As long as the null space 
of this linear differential expression can be analyzed in detail, it becomes possi- 
ble to transfer solutions, in fact, entire classes of solutions (e.g., rational, soliton, 
algebro-geometric solutions, etc.) between these evolution equations. In concrete 
applications, however, it turns out to be simpler to rewrite a relationship between 



two evolution equations, such as (1.4) and (1.28), in a form analogous to (1.6) and 
(1.2E), rather than analyzing the nullspaces of (2<p ± d x ) and (— d x + (ip x /ip)) in 
detail. These strategies relating (hierarchies of) evolution equations and their mod- 
ified analogs is not at all restricted to the Burgers and heat equations and the KdV 
and mKdV hierarchies, respectively, but applies to a large number of evolution 
equations including the Boussinesq fl29| , and more generally, the Gelfand-Dickey 
hierarchy and its modified counterpart, the Drinfeld-Sokolov hierarchy H|, the 
Toda and Kac-van Moerbeke hierarchies [|| , J33| , the Kadomtsev-Petviashvili and 
modified Kadomtsev-Petviashvili hierarchies [271] , J32|, etc. 

For simplicity we restrict ourselves to classical solutions throughout this note. 
The case of distributional solutions of Burgers equation is considered, for instance, 
in 0. 

Throughout this note we abbreviate by C p ' q (Q x A), ft C K", A C K open, neN, 
p, q G No, the linear space of continuous functions f(x,t) with continuous partial 
derivatives with respect to x = (x\, . . . , x n ) up to order p and q partial derivatives 
with respect to t. C p ' q (fl x A;R") is then defined analogously for f(x,t) G M n . 
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2. The Miura transformation 

We turn to the precise formulation of the relations between the KdV and the 
mKdV equation and omit details of a purely calculational nature. 

Lemma 2.1. Let if) = tp(x, t) > be a positive function such that ip G C 4: °(R x R), 
tpt G C h0 (R x R). Define 4> = i/} x /ip. Then 4> G C 34 (R x R) and 



mKdV(</>) = 9 X - P(V±)i/>)\ , (2.1) 



w/iere 

P(V) = 2V8 X - V x (2.2) 

and 

V ± =<f> 2 ±0 x . (2.3) 

Proof. A straightforward calculation. □ 

The application to the KdV equation then reads as follows. 

Theorem 2.2. Let V = V(x,t) be a solution of the KdV equation, KdV(V^) = 0, 
with V G C 3 ' 1 (R x R), and let ip > be a positive function satisfying ip G C 2 '°(R x 
R), ip t G C 1 '°(R x R) and 

ipt = P{V)i>, ~i>*x + Vip = 0, (2.4) 



with P(V) given by (2.2). Define cf> — ipx/ip an d V = 4* §x- Then V = 4> 2 + <fi x 
and <j) satisfies (f> G C 4,i (R x R) and the mKdV equation, 

mKdV(0) = 0. (2.5) 

Moreover, V satisfies V G C 3,1 (R x R) and the KdV equation, 

KdV(F) = 0. (2.6) 



Proof. A computation based on Lemma 2.1. □ 



Originally, Theorem |2J was proved in (2J (see also |21[, |g| (2J, |2§, |gl)) 
using supersymmetric methods. The above arguments, following ]27| in the context 
of the (modified) Kadomtsev-Petviashvili equation, result in considerably shorter 



calculations. The "if part" in Theorem 2.2 also follows from prolongation methods 



developed in |4Sf ■ A different approach to Theorem 2.2, assuming rapidly decreasing 
solutions of the KdV equation, can be found in Sect. 38 of §. 

Remark 2.3. The chain of transformations 

V ^ (j) ^ -0 V (2.7) 

reveals a Backlund transformation between the KdV and mKdV equations (V — > <fi) 
and two auto-Backlund transformations for the KdV (V — > V) and mKdV equations 
((f) — ♦ —0), respectively. 

Remark 2.4. For simplicity we assumed ip(x,t) > for all (x,t) G R 2 in Theo- 
rem [2.2| . However, as proven by Lax one can show that ip(x, to) > for some 
to 6 1 and all x £ R actually implies ip(x,t) > for all (x, t) G R 2 (see also 
|24[ ). Moreover, in case V(a;, to) is real-valued, we note that L(to)ip{x, to) = has 
a positive solution i/)(x, to) > if and only if the Schrodinger differential expression 
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L{to) = ~d 2 r + V(x, to) is nonoscillatory at ±00 (cf. pH ). While the system of 
equations (2.4) always has a solution ip(x,t), (cf. Lemma 3 in it is the addi- 

tional requirement ^(x, i) > for all (a;, t) £ R 2 w hich renders (and hence V) 
nonsingular. Without the condition ip > 0, Theorem 2.2 describes (auto)Backlund 
transformations for the KdV and mKdV equations with characteristic singularities 
(cf. H). 

3. The Cole-Hopf transformation 

Finally we return to relations ( |l.28[) , (|l.29|) , and Jl.3C|) . Since they are all proved 
by explicit calculations we may omit these details and focus on a precise formulation 
of the results instead. 

Lemma 3.1. Let ip — ip(x, t) > be a positive function with ip £ C 3 '°(R x (0, 00)), 
ipt £ C x,0 (Rx(0,oo)). Define V = -2vip x /ip withv > 0. Then V £ C* 2 < x (Rx (0, 00)) 
and 

—d + — 

ip x ip 2 



V t + VV X - vV xx = 2v ( --d x + ^)(ipt- vip xx ) (3.1) 



= -2vd x (^{ip t -vip xx ) S J. (3.2) 

The extension to the case with a potential term F in the heat equation reads as 
follows. 

Lemma 3.2. Let F £ C 1,0 (R x (0, 00)) and assume ip — ip(x,t) > to be a 
positive function such that ip £ C 3 '°(R x (0, 00)), ip t £ C 1,0 (R x (0, 00)). Define 
V = -2isip x /ip with v>Q. Then V G C 2 ' X (M x (0, 00)) and 

V t + VV X - vV xx + 2vF x = -2vd x {h^ - vip xx - Fip)^J . (3.3) 

We can exploit these relations as follows. 

Theorem 3.3. Let F G C X '°(R x (0, 00)) and v > 0. 
(i) Suppose V satisfies V G C 2 ' X (R x (0, 00)), and 

V t + VV X - vV xx + 2vF x = (3.4) 

for some v > 0. Define 



ip(x,t) = exp (—^j 



dyV(y,t)\. (3.5) 



Then ip satisfies < ip £ C 3 ' 1 ^ x (0, 00)) and 



-7 - i^»x* - FV») = C(t) (3.6) 



/or some x-independent C £ 
(nj Let ip > be a positive function satisfying ip £ C 3,0 (R x (0, 00)), ip t £ C 1,0 (R x 
(0, 00)) and suppose 

ipt = vip xx + Fip (3.7) 

for some v > 0. Define 

V = _2v^. (3.8) 



Then V £ C* 2 ' X (R x (0,oo)) safc^es (3.4). 
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Remark 3.4. One can "scale away" C(t) in Theorem |3.3| (i) by introducing a new 
function ijj. In fact, the function ip(x : t) = ip{x,t) exp(— ds C(s)) satisfies 

i> = H'xx + Fiji. (3.9) 

Remark 3.5. Using the standard representation of solutions of the heat equation 
initial value problem, 

ipt = vrfxx, i>(x,0) = ^o(x), (3.10) 

assuming 

ip eC(R), VoC^) < C*iexp(C 2 |x| 1+7 ) for |x| > i? (3.11) 

for some R > 0, Cj > 0, j = 1,2, and < 7 < 1, given by (cf. Q, Ch. 3; Q, 
Ch. V) 

^ *) = „, \u9 f d V exp{-(x - y) 2 /4vt)Mv) > 0, (3.12) 



the corresponding initial value problem for the Burgers equation 

V t + VV x -vV xx = 0, V(x,Q) = V {x), (3.13) 

reads 

v( _ / R dy (x - y)t^ exp ( - (2^ g d V V Q (r,) (x yfjAvt)^) 
{X ' ' J R dye^(-(2^^dr,V ( V )-(x-yn4ut)-^ > { ' > 

assuming Vq £ C(R) and 



V Q {x) > or 



dyV (y) 



= 0(M 1+7 ) for some 7 < 1. (3.15) 

\x\ — >oc 



Without going into further details we mention the existence of a hierarchy of 
Burgers equations which are related to the linear partial differential equations 

^ = <9™+V, neN (3.16) 

by the Cole-Hopf transformation V — —2tp x /tp (see ||). 

The multi-dimensional extension of Lemma |3.2| reads as follows. 

Lemma 3.6. Let F e C 1,0 (E n x (0, 00)) and assume ip — ip(x,t) > to be a 
positive function such that ip G C 3 < (R" x (0,oo)), ip t £ C 1 '°(R n x (0, 00)). Define 
V = -(2i//a)Vhi(i/i) with a G K\{0}, v > 0. Then V G C 24 (K" x (0,oo);R Tl ) and 

V t + a(V- X/)V-vAV+— VF= -— V ( - U t - vM> - Ftp) ) . (3.17) 

a a \ip J 

Our final result shows how to transfer solutions between the multi-dimensional 
Burgers equation and the heat equation. 

Theorem 3.7. Let F £ C lfi (R n x (0,oo)), a G R\{0}, and v > 0. 
(%) issume i/iat V G C 2 - 1 ^™ x (0,oo);R n ) satisfies 

V = V$ (3.18) 

/or some potential <£> £ C 3,1 (R™ x (0, 00)) and 

V t + a(V ■ V)V - vAV + — VF = 0. (3.19) 

a 
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Define 

Then tp S C 3,1 (R n x (0,oo)) and 



h^t-v^-F^) = <7(t), (3.21) 
for some x -independent C £ C((0, oo)). 

(ii) Let xjj > be a positive function satisfying ip G C 3 '°(R™ x (0, oo)), ipt £ 
C 1,0 (M. n x (0, oo)) and suppose 

vb t = v^ + Fvb. (3.22) 



Define 



V = -— Vln(^). (3.23) 
a 



Then V e C 24 (M n x (0,oo);R") satisfies (3.19) 



Acknowledgments. We thank Mehmet Unal and Karl Unterkofler for discus- 
sions on the Burgers and (m)KdV equations, respectively. 
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